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nificant, short-timescale deviation from the power-law flux decline expected in 



The optical afterglow of Gamma-Ray Burst (GRB) 000301C exhibited a sig- 



the standard synchrotron shock model. Garnavich, Loeb & Stanek found that 
this deviation was well-fit by an ad hoc model in which a thin ring of emission 
is microlensed by an intervening star. We revisit the microlensing interpreta- 
tion of this variability, first by testing whether microlensing of afterglow images 
with realistic surface brightness profiles (SBPs) can fit the data, and second by 
directly inverting the observed light curve to obtain a non-parametric measure- 
ment of the SBR We find that microlensing of realistic SBPs can reproduce the 
observed deviation, provided that the optical emission arises from frequencies 
above the cooling break. Conversely, if the variability is indeed caused by mi- 
crolensing, the SBP must be significantly limb-brightened. Specifically, > 60% 
of the flux must originate from the outer 25% of the area of the afterglow image. 
The latter requirement is satisfied by the best fit theoretical SBP. The underlying 
optical/infrared afterglow lightcurve is consistent with a model in which a jet is 
propagating into a uniform medium with the cooling break frequency below the 
optical band. 
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Introduction 



The afterglows of gamma-ray bursts (GRBs) are observed in the X-ray, optical, near- 
infrared, and radio, and appear to be well-described by the synchrotron blast-wave model in 
which the source ejects material with a relativistic bulk Lorentz factor, driving a relativistic 
shock into the external medium (see Van Paradijs et al. 2000; Piran 2000 and references 
therein). There is mounting evidence from the observed steepening of afterglow light curves 
that these ejecta are in many cases mildly to highly collimated, with opening angles ~ 3-30° 
(Harrison et al. 1999; Frail et al. 2000, 2001). Global fitting of the afterglow light curves 
over many decades in time and frequency, in the context of this model, can be used to derive 
constraints on the physical parameters of the model, i.e., the energy and opening angle of 
the jet, the external density, the magnetic field strength and the energy distribution of the 
electrons behind the shock (Wijers & Galama 1999; Freedman & Waxman 2001; Panaitescu 
& Kumar 2001). 

In this model, the image of the afterglow is expected to appear highly limb brightened at 
frequencies above the peak synchrotron frequency z/ m , but more uniform at frequencies < u m , 
especially below the self-absorption frequency v a (Waxman 1997; Sari 1998; Panaitescu & 
Meszaros 1998; Granot, Piran, & Sari 1999a,b; Granot & Loeb 2001). A measurement of 
the surface brightness profile (SBP) at several frequencies would thus provide an important 
test of the model. For typical parameters, the afterglow image expands superluminally, and 
has an angular radius (9(/ias) a few days after the GRB. As pointed out by Loeb & Perna 
(1998), this is of the same order as the angular Einstein ring radius of a solar mass lens at 
cosmological distances, 



where M is the lens mass, and D = D 0S D i/Di s , and D os , D Q \ and D\ s are the angular 
diameter distances between the observer-source, observer-lens, and lens-source, respectively. 
Thus lensing by a star along the line of sight will produce a detectable magnification over 
the course of a few days to weeks. The probability that any given GRB will be microlensed 
is ~ l%b 2 , where b is the angular separation between the GRB and lens in units of $e- Since 
the original suggestion by Loeb & Perna (1998), this application of microlensing has been 
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studied by several authors (Mao & Loeb 2001; Koopmans & Wambsganss 2001; Granot & 
Loeb 2001; Gaudi & Loeb 2001; Ioka & Nakamura 2001). 

The well-sampled optical afterglow of GRB 000301C exhibited a short time scale, nearly- 
achromatic deviation from the nominal power-law flux decline. Garnavich, Loeb & Stanek 
(2000; hereafter GLS) found that this deviation is well-fit by a model in which the afterglow 
was microlensed by an intervening star. However, the model they adopted for the SBP was 
somewhat crude: they assumed the emission arose solely from an outer ring of fractional 
width W, with W ~ 10% providing the best fit to the data. Recently, Panaitescu (2001) 
argued that, if one adopts realistic SBPs, microlensing cannot explain the data. Here we 
revisit the microlensing interpretation of this deviation. We fit the optical/infrared dataset 
of GRB 000301C to a double-power law plus microlensing model, adopting realistic SBPs, 
and, under the assumption that the observed deviation is indeed due to microlensing, we 
directly invert the light curve, obtaining a non-parametric measurement of the SBP, and 
compare this measurement to theoretical expectations. 

2. Data 

We will be considering the same dataset as GLS, namely 104 photometric data points 
in the UBVRIJK-b&nds distributed as (6,18,8,46,16,3,7). These data, shown in Figure 
1, are taken from the compilation of Sagar et al. (2000) with additional photometry from 
Stanek et al. (2000). We will not consider the data in the radio regime (Berger et al. 2000), 
as the majority is not contemporaneous with the optical deviation. For more details on the 
GRB itself and the optical, infrared and radio data, see GLS, Masetti et al. (2000), Sagar 
et al. (2000), Berger et al. (2000) and Rhoads & Fruchter (2000). 

3. Theoretical Surface Brightness Profiles 

The magnitude, shape, and duration of the microlensing signal depends sensitively on 
the SBP of the afterglow (Loeb & Perna 1998; Mao & Loeb 2001; Granot & Loeb 2001). The 
SBP, in turn, depends on the observed frequency and the physical parameters of the shock. 
Thus detailed measurements of a microlensed afterglow light curve can be used to constrain 
the physical parameters of the afterglow. However, the expected SBPs in the relativistic 
blast-wave model can also be calculated a priori under various assumptions (Waxman 1997; 
Sari 1998; Granot, Piran & Sari 1999a,b; Granot & Loeb 2001; Ioka & Nakamura 2001; 
Panaitescu 2001). Thus one can alternatively use the theoretical SBPs to calculate the 
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Fig. 1.— UBVRIJK photometry of GRB 000301C as a function of days from the GRB. 
Points have been offset by the indicated amount for clarity. The dashed line is the best-fit 
double power law (with no lensing), while the solid line is the overall best-fit microlensing 
model, where the SBP has been determined from direct inversion. The inset shows the 
i?-band data only. The dotted line is the best-fit microlensing model with theoretically 
calculated SBP, in this case for k = and v > v c . 
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expected microlensing signature, thereby confirming or refuting the microlensing interpre- 
tation 1 . This was the approach taken by Panaitescu (2001). We adopt both approaches 
here. 

We calculate SBPs using the method described in Granot & Sari (2001) and Granot & 
Loeb (2001). Briefly, the hydrodynamics are described by the Blandford-McKee (1976; BM 
hereafter) self-similar solution, assuming a power-law external density profile, p(r) oc r~ k , 
and a power-law number versus energy distribution of electrons with index p, dN e /d , ~f oc 7 _p , 
just behind the shock, which thereafter evolves due to radiative and adiabatic losses. We 
integrate the emission over the entire volume behind the shock front. A few days after the 
burst the afterglow is typically in the slow cooling regime (i.e. a typical electron cools on 
a timescale larger than the dynamical time) and the ordering of the break frequencies is 
v a <v m < v ci where v c is the cooling frequency (Sari, Piran, & Narayan 1998; Granot, Piran 
& Sari 1999a; Chevalier & Li 2000). We assume that this is the case for GRB 000301C, 
as indicated by Berger et al. (2000). We consider only the SBPs from frequencies v > u a , 
since after a few days the optical band is typically well above u a (and usually also above 
v m ). The spectrum consists of several power-law segments (PLSs), where F v oc v"t a (see 
Figure 1 in Granot & Loeb 2001). Due to the self-similar nature of the hydrodynamics, the 
normalized SBP within a given PLS is independent of time (though it changes considerably 
between different PLSs). Thus, for each k, there are 3 relevant forms for the SBP, each 
corresponding to a different PLS: f3 — 1/3 (u a < v < u m ), ft — (1 — p)/2 [y m < v < u c ) and 
(3 = —p/2 [y > u c ). The corresponding values of the temporal index a for k = (k = 2) are 
a = 1/2 (a = 0) for (3 = 1/3, a = -3(p - l)/4 (a = -(3p - l)/4) for f3 = (1 - p)/2 and 
a = -(3p- 2)/4 (for both k = and k = 2) for (3 = -p/2. 

The expected SBPs are shown in Figure 2, for a uniform external medium k — 0, and 
a stellar wind environment k = 2. Note that the SBP becomes more limb-brightened as the 
frequency goes to a higher PLS, i.e. the spectral slope f3 is reduced, or as the value of k is 
reduced. 



1 Since both the unlcnscd light curves and the magnification history depend on the physical assumptions 
in the afterglow model, a poor fit for a certain microlensing model might suggest an unrealistic afterglow 
model, rather than refute the microlensing interpretation altogether. 
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Fig. 2. — SBPs as a function of normalized radius. The points are the SBPs determined from 
direct inversion, with la errors (solid) and 3a errors (dotted). The curves are theoretically 
calculated SBPs for various frequency regimes (see Granot & Loeb 2001). (a) Uniform 
external medium, k = 0. (b) Stellar wind environment, k = 2. 
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4. Analysis and Model Fitting 



We fit the data of the flux as a function of time to the form of two asymptotic power-laws 
which join smoothly: 



where % is the "break time," Fq^ is the unlensed flux at the point where the asymptotic 
power laws before and after the break meet, ol\ and a<i are the indices of the flux decline, and 
s is a sharpness parameter; the larger the value of s, the sharper the break. The magnification 
[i due to microlensing depends on the impact parameter b, the radial SBP of the afterglow 
image, and the angular size of the afterglow in units of Oe, 



Here 5 = (5 — fc)/2(4 — k), tday is the time from the GRB in days and R Q is the value of 
R s after 1 day. The light curve exhibits an achromatic break at t h ~ 7 days. Berger et 
al. (2000) find that this break can be well explained by jet with k — 0, while Panaitescu 
(2001) argues that a jet alone cannot produce the required steepening of the light curve, and 
invokes a double power law electron distribution, which causes additional (though chromatic) 
steepening of the light curve if the synchrotron frequency corresponding to the break in the 
electron distribution crosses the optical band around t\>. Either way, we would expect the 
SBP at t~ tb to be different than at t~%. However, our best fit models give % ~ 8 days, while 
the majority of the constraints on the SBP come from times t~8 days, where the spherical 
self-similar BM solution should still be applicable, resulting in the self-similar SBPs described 
in §3. We will therefore adopt the simpler (and self consistent) working assumption that 
the SBP, when normalized to its average value, as a function of radius normalized to R s , is 
independent of time within a given PLS. 

The model in equation (2) is a function of jV par = A + N t + 2 free parameters: four 
parameters describe the shape of the broken power-law (cci, a 2 , s, ty,); one parameter pro- 
vides the normalization Fq jV for each of iVj frequency bands 2 and two additional parameters 



2 Standard afterglow models predict Fo iV oc . However, extinction from the Milky Way and possibly 
the host galaxy itself will induce an (uncertain) curvature in the spectrum. Indeed, as noted by Rhoads 
& Fruchtcr (2000); Jensen et al. (2001) and Panaitescu (2001), the early optical/infrared data for GRB 
000301C show evidence for a curved spectrum, which may be due to extinction or is perhaps intrinsic. For 
simplicity, we will use a free parameter for each Fo t „. We have also performed tests where we fit the data 
assuming the form F QjV oc v 13 . In general we find worse fits, though our conclusions regarding the SBPs are 
essentially unchanged. As we show, the results for the i?-band data are the same as for the full dataset, 
indicating that our assumption of independent F 0j „ does not bias our results. 
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describe the microlensing effect (b, R ) for a given SBP. Considering all the data for GRB 
000301C in Figure 1, there are 11 parameters without microlensing, and 13 free parameters 
with microlensing using theoretically-calculated SBPs. For comparison purposes, we also fit 
the data to the model adopted by GLS: we assume that the source emission is confined to a 
thin ring of fractional width W. However, we extend the model of GLS by allowing the image 
interior to the ring to have a relative surface brightness C referenced to the outer ring (C = 
reduces to the GLS model). This model introduces two additional free parameters. Finally, 
we consider a non-parametric fit ("direct inversion") to the SBP. We divide the image into 
Nun equal-area annuli, and find the fraction fa of the total flux contributed by each bin i. In 
this case the magnification is given by /i — J2i fiVu where ^ is the magnification of annulus 
%. This adds an additional N^ n — 1 parameters (due to the constraint that £V fi — 1). We 
find the best-fit solution by minimizing x 2 with respect to all of these parameters. We define 
the 1(7 errors on these parameters as the projection of the A% 2 = 1 hypersurface on the 
parameter axes, where 

A X 2 = ^Pjff , (4) 

Xmin i s the minimum x 2 f° r a given model, and dof = 104 — N par is the number of degrees 
of freedom for the 104 data points. We normalize Ax 2 by the factor Xmin/dof because we 
believe that the errors on the data points are likely underestimated, resulting in inflated 
values of x 2 ■ Errors on the fit parameters determined using these inflated values of x 2 would 
be significantly underestimated. 



5. Results 

Figure 1 shows the best fit for the double power-law model with no lensing. The best- 
fit parameters and la errors for (a^, «2, s, %) are tabulated in Table 1. The fit is poor: 
X 2 = 240.7 for 93 dof. Residuals from the broken power-law model are shown in Figure 3; 
the systematic deviations from this model are apparent. 

When microlensing is included, the fit improves considerably. For the extended GLS 
model, we find that the data is best explained by emission solely from an outer ring of 
relatively small fractional width. Specifically, we find best-fit parameters, W = 0.13^oo7> 
C = 0+ 02 for k = and W = 0.16lg;g|, C = 0+ 03 for k — 2. The remaining parameters 
and the best-fit x 2 are shown in Table 1. The parameters we find differ from those of GLS, 
and we find a x 2 /dof smaller than they report. However, when we evaluate their model fit, 
we recover their x 2 /dof, indicating that their fit was likely a local, not global, minimum. For 
the fits with microlensing, we find that the sharpness parameter s is not well constrained, 
but is always s~5, indicating a sharp break. To avoid excessive covariances with other 
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parameters, we fix s = 20 for the remaining fits, thus reducing N pai by one. The results 
are essentially identical for other values of s~5. For a uniform external medium (k = 0), 
for which one generally predicts more limb-brightened SBPs, we find \ 2 = 190.0, 166.1, and 
154.7 for v a < v < u m , v m < v < u c , and v > v c , respectively. For the case of a wind-like 
medium (k = 2), the fits are generally worse: \ 2 = 216.8,197.9, and 155.5 for the same 
frequency ranges. The deviations of these fits from the double power-law model are shown 
in Figure 3. It is clear that a significantly limb-brightened profile is required to fit the sharp, 
large deviation at t ~ 4 days. For v > u c , both wind and uniform external media provide 
satisfactory fits to the data. For u, m < u < v c , the uniform medium provides a marginal 
fit. Frequencies below the peak synchrotron frequency u m , or u m < v < u c with a wind-like 
medium (k — 2), result in SBPs that are too uniform to fit the data. 

For the non-parametric SBP fit to the data, we impose the condition that fa > 0, and 
we increase the number of bins A^ bin in the recovered SBP until x 2 changes by < 1. We 
find that N hin = 4 is optimal. We find \ 2 = 148.5 for k = 0, and \ 2 = 146.9 for k = 2, 
for 89 dof. The inferred surface brightnesses are shown in Figure 2, the fit for k = is 
shown in Figure 1, and Figure 3 shows the deviations of both the k = and k = 2 models 
from the best-fit double power law. The best-fit values and la errors for the parameters a±, 
a 2 , tb, Ro, and b are given in Table 1 for all viable (Ax 2 < 20) fits. Although the values 
of x 2 /d°f for the direct inversion formally indicate poor fits, it is clear from Figures 1 and 
3 that microlensing reproduces the observed light curve shape quite faithfully, at least for 
limb-brightened sources. It is likely that underestimated errors contribute significantly to 
the inflated x 2 - 

We have also fit the i?-band data only using both the SBP determined from direct 
inversion and the theoretically-calculated SBPs. The SBP inferred from direct inversion is 
consistent with that found from fitting all the data (with larger errors, of course). The 
fit is better: x 2 /dof = 1.23 for 37 dof. Fits using the theoretical SBP yield the same 
conclusions as the fits to the entire data set: the SBPs for v > v c provide fits that are 
acceptable (x 2 /dof = 1.24 and 1.26 for k = and k = 2 respectively, for 40 dof). The fit for 
v m < v < u c , and k = is marginal: xVdof = 1.68. All other fits are considerably worse: 

x 2 /dof > 2. 



6. Discussion and Conclusions 

We find that the early-time (t < tb) flux decline is well-constrained, ol\ = —1.05 ± 0.02. 
The slope of the optical spectrum at this time is f3 — —0.9, taking into account Milky Way 
(MW) extinction only, but can be as shallow as [3 = —0.7 if SMC-type host extinction is 
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Days Since GRB 

Fig. 3. — In each panel, the points show the binned residuals from the best-fit double power- 
law fit (without lensing) in magnitudes. The curves show the difference between microlensing 
models with various SBPs and the double power law. (a) Theoretically calculated SBPs for 
k=0. (b) Theoretically calculated SBPs for k=2. (c) SBPs determined from direct inversion. 



- 11 - 



Model 


— Cti 


— a 2 


s 




6 


Rq 


X 2 /dof 


No /i-lensing 

k = 0, u m < v < u c 
k = 0, v > v c 
k = 2, v > u c 
k = 0, DI 2 
k = 2, DI 2 
k = 0, GLS 3 
fc = 2, GLS 3 


_ — in no 
± - u ' -0.02 

l 06 +0 - 02 

± - uu -0.02 
i - u ^-0.02 

1 06 +a02 

1 - uu -0.02 
-■-• uo -0.02 

1 06 +a02 

i - uu -0.02 
i - uo -0.02 


9 79+0.11 
z - ' z -0.10 

901 + U.ZO 
z - oi -0.10 

2 62 +(U0 

z - uz -0.11 

2 63+ 011 

z - uo -0.11 

2 60 +an 

z - uu -o.io 

9 ci +0.11 
Z - Ji -0.10 
9 71 +0.10 
z - ' x -0.09 
9 7 r+0.10 
z - ' °-0.10 


14+1.37 

-J-^-o.gi 
20 1 

20 1 

20 1 

20 1 

20 1 

iJ - J -10.0 

13.81& 


4.8olg : S 

1 4+U.54 
7.89lgjg 

s.oi±8iJ 

7 OO+0.45 
' -°°-0.41 

8.361^ 
7.68lgjJ 


u - z ^-0.24 

0.651™ 
0.601™ 
62 +0 ' 24 

u - uz -0.19 

0.4418^ 

o.8718;1S 
o.92i8:2 


_ _ _ _i_n 1 1 

n 1 9+0- 11 

u - iz -0.12 

2Q+ - 06 

u -^ y -0.06 

n 23-°- 05 

u - zo -0.05 
u - zy -0.06 

n i8+ - 13 

n 41 +0.07 
u -^ i -0.08 
n 07+O.O6 
u -° ' -0.07 


240.7/93 
166.1/92 
154.7/92 
155.5/92 
148.5/89 
146.9/89 
147.3/89 
146.9/89 



Table 1 Fit parameters and la errors for GRB 000301C. 



1 Fixed. 

2 Fit with the SBP determined from direction inversion. 

3 Fit using a modified version of the SBP model adopted by GLS. 

included (Rhoads & Fruchter 2000; Jensen et al. 2001). For v > u c (the SBPs most consistent 
with the data), the value of a\ implies p ~ 2.07 and thus (3 = —1.03, which is only slightly 
steeper than the observed spectrum for MW extinction. For v m < v < v c and k — 0, 
cti = —1.05 implies p = 2.4 and thus (3 = —0.7, consistent with observations z/SMC-type 
host extinction is adopted. If the break at t ~ 8 days is due to a jet, then simple analytic 
models predict that the change in temporal flux index a should be Aa = (p + 3)/4 for k = 
and v m < v < v c or Act = (p + 2)/4 for v > v c and k = or 2 (Sari, Piran, & Halpern 
1999; Panaitescu & Kumar 2001). Kumar & Panaitescu (2000) find, using a semi-analytic 
model, that the jet break is rather smooth, especially for a wind environment (k = 2), 
and therefore inconsistent with the sharp break that we infer. However, initial results of 
numerical calculations of the jet break (Granot et al. 2001), based on a 2D hydrodynamic 
simulation of the jet for k — 0, indicate that the break can be rather sharp (with s ~ 4.5 
for an observer along the jet axis), and that a 2 is smaller by ~ 0.35 compared to the simple 
analytic predictions mentioned above, making Aa larger by a similar factor. Nevertheless, 
the sharpness of the break still presents a serious problem for the model where the jet 
propagates into a stellar wind (k = 2). Also, the SBP is expected to become more uniform 
at t > % compares to t < % (Ioka & Nakamura 2001; Panaitescu 2001). Since we assume 
that the form of the SBP relevant for t < % holds all along, we underestimate the value of 
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«2, and overestimate the value of Act accordingly 3 . This effect is expected to be larger for 
v > v c , where the SBP is more limb brightened to begin with. The values we obtain for Ao 
may therefore serve as upper limits (especially for v > v c ). For the SBP fit for k — and 
u m < v < u c , we find Act < l-2lg 1 whereas based on the electron index derived from ai, the 
prediction is Act ~ 1.7. For v > v c and k = 0, we find Act < 1.56 ±0.11, which is reasonably 
consistent with the prediction Aa pa 1.37. 

Thus, assuming the break in the light curve is due to jet effects, the only scenario which 
provides a satisfactory explanation of the data is v > v c and k — 0. For v < u c , the inferred 
SBP is only marginally consistent with the expected profile, and the observed change in the 
flux indices is smaller than the theoretical expectation. For v > v c and k = 2, the observed 
break is too sharp compared to theoretical predictions. We note that our best fit model 
parameters, v > v c and k — 0, are also favored by Panaitescu (2001), and said to provide a 
good fit to the data by Berger et al. (2000). 

Panaitescu (2001) argued that the expected SBPs are too uniform to reproduce a signif- 
icant microlensing deviation. This conclusion appears to be partially due to the fact that he 
found a best-fit impact parameter of b ~ 2, whereas we find b = 0.44-0.65 (and b = 0.65 for 
our best fit model). All else equal, more uniform SBPs will result in microlensing deviations 
that are shallower and broader than more limb-brightened SBPs. Thus, to reproduce the 
observed deviation, a more uniform source must have a smaller impact parameter, a trend 
we find empirically when fitting the theoretical SBPs. Thus, it is unclear why Panaitescu 
(2001), who argued that the expected SBPs should be more uniform than that found by 
GLS, derived a value of b that is twice as large (b ~ 2 versus b ~ 1). 

We conclude that the microlensing of realistic SBPs can explain the optical/infrared 
light curve of GRB 000301C. Fitting the data with a non-parametric SBP reveals that the 
afterglow image must be significantly limb-brightened. Specifically, ~ 60% of the flux must 
originate from the outer 25% of the area of the image (with a significance of 3a). The 
recovered SBP is marginally consistent (at the ~ 3<r level) with the profile expected for a 
uniform external density profile (k = 0) and frequencies v m < v < u c . It is consistent at the 
~ 2a level with the emission expected from frequencies u > u c and both uniform (k = 0) 
and wind-like (k = 2) external media. 



3 This effect may be seen in Table 1, where 012 is larger for v m < v < v c (for which the SBP is more 
homogeneous), compared to v > v c , for k = or 2, (for which the SBP is more limb brightened). 
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